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ABSTRACT 

We investigate the utility of the Tunable Filters (TFs) for obtaining flux calibrated emission line maps 
of extended objects such as galactic nebulae and nearby galaxies, using the Optical System for Imaging 
and low Resolution Integrated Spectroscopy (OSIRIS) at the 10.4-m Gran Telescopio Canarias (GTC). 
Despite the relatively large field of view of OSIRIS (8' x 8'), the change in wavelength across the field 
(~ 80 A) and the long-tail of the TF spectral response function are hindrances for obtaining accurate 
flux calibrated emission-line maps of extended sources. The purpose of this article is to demonstrate 
that emission-line maps useful for diagnostics of nebula can be generated over the entire field of view of 
OSIRIS, if we make use of theoretically well-understood characteristics of TFs. We have successfully 
generated the flux-calibrated images of the nearby, large late-type spiral galaxy MlOl in the emission 
lines of Ha, [Nii]A6583, [Sii]A6716 and [Sii]A6731. We find that the present uncertainty in setting 
the central wavelength of TFs (^lA), is the biggest source of error in the emission-line fluxes. By 
comparing the Ha fluxes of Hii regions in our images with the fluxes derived from Ha images obtained 
using narrow-band filters, we estimate an error of ~ 11% in our fluxes. The flux calibration of the 
images was carried out by fitting the Sloan Digital Sky Survey (SDSS) griz magnitudes of in-frame 
stars with the stellar spectra from the SDSS spectral database. This method resulted in an accuracy of 
3% in flux calibration of any narrow-band image, which is as good as, if not better than, that feasible 
using the observations of spectrophotomctric standard stars. Thus time-consuming calibration images 
need not be taken. A user-friendly script under the IRAF environment was developed and is available 
on request. 

Subject headings: galaxies: photometry — ISM: HII regions — methods: data analysis — techniques: 
image processing 



1. INTRODUCTION 

Optical System for Imaging and low Resolution Inte- 
grated Spectroscopy (OSIRIS) is an imager and spec- 
trograph for the optical wavelength range, located in 
the Nasmyth-B focus of the 10 . 4-m Gran Telesco - 
pio Ganarias (GTC) (|Cepa et all [20051 : ICepal [20ll . 
Apart from the standard broad-band imaging and long- 
slit spectroscopy, it provides additional capability such 
as narrow-band tunable filter (TF) imaging, charge- 
shuffling and multi-object spectroscopy. OSIRIS covers 
the wavelength range from 0.365 to 1.05 /zm with an un- 
vignetted field of view (FoV) of 7.8 x 7.8arcmin^, for 
direct imaging. Narrow-band imaging is made possible 
through the use of a TF, which is in essence a low reso- 
lution Fabry-Perot etalon. The filter is tuned by setting 
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a specific separation between the optical plates of the 
Fabry-Perot. OSIRIS has two different TFs: one for the 
blue range (3750-6750 A; yet to be commissioned), and 
another for the red range (6510-9350 A). In this work, 
we describe the results for the latter (red) etalon. 

The relatively large FoV on a 10-m telescope, com- 
bined with the tunable nature of the filters, opens up 
a new method of observation, which can be useful in a 
variety of astrophysical contexts, both galactic and extra- 
galactic. The use of narrow-band filters, adequately de- 
fined to map certain spectral regions, allows the detailed 
study at seeing-limited resolution, of the spatial distri- 
bution of the ionized gas properties and recent star for- 
mation, and of the associated stellar populations. Flux 
calibrated images in the bright lines of hydrogen, oxygen, 
nitrogen and sulphur in the optical, allow us to map, 
among other parameters, the metallicity gradients in 
galaxies using imaging techniques rather than the time- 
consur ning spectroscopic technique s that are in use cur- 
rently (jRosales-Ortega et al.|[2010r ). This has motivated 
recently built large telesc opes such as the 11-m South 
African Large Telescope (iRangwala et al.l I2008D or the 
Magellan Baade 6.5-m telescope ( Veilleux et al .1120101 ) to 
have TFs for imaging. We are presently engaged in a 
survey of galaxies in the local universe (LUS^= Local 
Universe Survey) that includes all galaxies inside a vol- 
ume of 3.5 Mpc in radius that are visible from La Palma. 

"^ http://www.inaoep.mx/~gtc-lus/ 
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IFU observations of these galaxies present two problems: 
IFU observations would represent a huge drop in spatial 
resolution and, with a median D25 of 4.5', they are too 
la rge to be observed even w ith the largest available IFU. 

iLara-Lopez et al.l (|201Clf ) used template spectra typ- 
ical of star-forming galaxies to estimate the errors in 
derived flux for extragalactic point-like sources on the 
OSIRIS TF images. However, photometric accuracy 
of the TF images of extended objects using observed 
data is yet to be evaluated. In this study, we explore 
the capabilities of TF for obtaining flux calibrated im- 
ages using observations acquired with OSIRIS TF. The 
OSIRIS etalon is placed in a converging beam, which 
results in the variation of th e effective wavelength over 
the FoV ()Bornfc WolJll980D . For the red etalon, the 
change is as much as 80 A over the 4' radius, thus 
limiting the monochromati c FoV to less than 2' radius 
()Mendez-Abreu et al.|[201ll ). In addition, the TF spec- 
tral response function (see Eq. [2]) has a long tail, which 
difficults the flux calibration of a faint line in the neigh- 
borhood of a bright one. For example Ha line contam- 
inates even if the TF is tuned to detect the [Nii]A6583 
line. Given that these characteristics of TF are well- 
understood, flux calibrated monochromatic images over 
the entire FoV can be obtained using a tailor-made post- 
observation reconstruction software. 

The fact that the filters are tunable, in addition, cre- 
ates its own problems for absolute flux calibration due to 
the non-standard nature of the filters. Observing spec- 
trophotometric standard stars at each tuned wavelength 
is a solution, but it is costly in terms of telescope time. 
Alternatively, small telescopes could be used to carry 
out spectrophotometric observations of the in- frame field 
stars. The recent Sloan Digital Sky Su rvey (SPSS) spec- 
trophotometric survey of field stars ([Beers et al.l 120061 : 
lAbazajian et al. 2009) offers an attractive solution to the 
problem, that does not require any telescope time for cal- 
ibration purposes. In this paper, we explore the use of 
this database to spectrophotometrically calibrate the in- 
frame SDSS photometric stars. 

In §2, we describe the data used in this work. The 
method we followed for reconstructing a monochromatic 
image and its implementation are described in §§3 and 4, 
respectively. The sensitivity of the reconstructed emis- 
sion line fluxes to various observational parameters is 
studied using simulated images in §5 where we also sug- 
gest guidelines for planning imaging observations with 
OSIRIS/TF. The flux calibration technique and accuracy 
are described in §6. Ha fluxes of selected H 11 regions in 
the reconstructed image are compared to an Ha image 
obtained in a traditional way in §7. The results of this 
study are summarized in §8. 



2. DESCRIPTION OF THE DATA FOR RECONSTRUCTION 

2.1. The OSIRIS detector and image formats 

OSIRIS uses two CCDs of 2048 x 4102 pixel format 
to cover its total FoV of 8' x 8', with a physical gap 
of ~ 9" between the two CCDs. Each one has 50 ad- 
ditional overscan pixels at the beginning of the array, 
thus resulting in a format of 2098 x 4102 for each im- 
age. The optica l center of OSIRIS lie s in th e central gap 
as described in iMendez-Abreu et al.l ()2011l) and also in 







Table 


1 




GTC 


OBSERVING BLOCKS SUMMARY 


Observing Block Ac 


FWHM 


Order Sorter 


Exp. Time 




(A) 


(A) 


Filters 


(seconds) 


Ha-(-[Nll] 


6528 


18 


f648/28 


2 X 180 




6548 


18 


f648/28 


2 X 180 




6568 


18 


f648/28 


2 X 180 




6588 


18 


f657/35 


2 X 180 




6608 


18 


f657/35 


2 X 180 




6628 


18 


f657/35 


2 X 180 




6648 


18 


f657/35 


2 X 180 




6668 


18 


f657/35 


2 X 180 




6688 


18 


f666/36 


2 X 180 


[Sii] 


6696 


18 


f666/36 


3 X 180 




6716 


18 


f666/36 


3 X 180 




6736 


18 


f666/36 


3 X 180 




6756 


18 


£666/36 


3 X 180 




6776 


18 


£680/43 


3 X 180 




6796 


18 


£680/43 


3 X 180 




6816 


18 


£680/43 


3 X 180 



the OSIRIS User Manual*. Astrometric calibration is 
required before being able to mosaic the images regis- 
tered in the two CCDs. By performing astrometry of 10 
stars on each CCD that are distributed over the entire 
FoV, we found that the image scale is slightly different 
for the two CCDs. The left and right CCDs have values 
of 0.127"/pixel and 0.129"/pixel, respectively, with an 
average value of 0.128"/pixel. This average value agrees 
well with the theoretically expected value of the plate 
scale for the GTC/OSIRIS instrument p arameters listed 
in Table 1 of iMendez-Abreu et al.l ()2011[ ). 

2.2. OSIRIS TF emission-line imaging of MlOl 

We illustrate the reconstruction technique using ob- 
servations of MlOl aimed at obtaining continuum-free 
images of the Ha, [N ii]A6583, [S ii]A6716 and [S ii]A6731 
emission lines over a FoV of 8' x 8' . The observations were 
carried out at the GTC in two observing runs on June 
22 and 26 in 2009, the former for the [S ii]AA6716,6731 
lines and the latter for the Ha, [Nii]A6583 lines. Ta- 
ble 1 describes the details of these two observing runs. 
The telescope pointing remained the same for all images 
(wavelengths) constituting a TF scan. The telescope po- 
sition was then dithered by 6" and the entire sequence 
was repeated. This was done in order to facilitate re- 
moval of any detector artefacts. For a given region in the 
galaxy, the two dithered TF images have slightly different 
wavelength and hence different response for the detection 
of a line. Thus it is important to take into account these 
response differences before coadding them. We handled 
the dithered image sets as independent sets of data and 
used them to estimate internal errors in the fiux cali- 
bration (see §6). For the Ha-|-[Nii] scan, observations 
were carried out at two dithered positions (PI and P2, 
henceforth), whereas for the [S 11] scan, three dithered 
positions were used (PI, P2 and P3, henceforth). 

Initial reduction of the images for correcting for BIAS 
was done independently for the two CCDs constituting 
an image. Reliable fiat-field frames were not available 

* http://www.gtc.iac.es/en/mcdia/documcntos/OSIRIS-USER- 
MANUAL.vl.l.pdf 
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for the set of data we were analyzing, hence no attempts 
were made to correct for this. Typical pixel-to-pixel er- 
rors in response are expected to be less than 1% for the 
CCD chips used. Stars in the field of each of the two 
CCDs were independently analyzed to obtain the astro- 
metric solution between the pixel and equatorial coordi- 
nates, using the SDSS coordinates of around 10 stars in 
the observed field as reference. The two CCDs were then 
stitched together to form a single image, free of all geo- 
metrical distortions. The procedure is repeated for the 
other scans. 

3. MONOCHROMATIC IMAGE RECONSTRUCTION TECHNIQUE 

3.1. Basic formulae for tunable filters 

The effective wavelength A^ at a distance r from the 
optical center of the tunable filter (TF) ch anges following 
the law (,Mcn dez-Abreu et al.|[20Tll : lBecke rs 1998): 



Ar - — 



A. 



Vl + 6.5247 X 10- 



9^-2 ' 



(1) 



where Ac is the wavelength at the optical center. The 
coefficient of r^ is the inverse square of the effective fo- 
cal length of the camera, with r expressed in physical 
units of pixels of 15/im in size. Due to the geometrical 
distortions of the images, a given Ar moves to a differ- 
ent radial distance r' in the astrometrically calibrated 
images. We used the astrometric solution to map the 
new wavelengths A,.' at each astrometrically calibrated 
image with resampled pixels of 0.125" in size. This cor- 
rection, which is not symmetric around the optical cen- 
ter, increases quadratically with radial distance from the 
optical center, reaching values of ~ 2 and '^ 5 A at 4' 
away from the optical center on the left and right CCDs, 
respectively. 

For a TF of Full Width at Half Maximum FWHM, the 
filter response for a certain wavelength A (e.g. the Ha 
line) in a pixel where the nominal wavelength is A^, is 
given by: 



Rx{r) 



2(A- Ar) 



FWHM 



(2) 



This equation is an approximation valid when the free 
spectral range (i.e. the wavelength separation between 
adjacent orders) is much larger than the spectral pu- 
rity (i.e. the smallest m easurable wavelength difference) 
— see, e.g.. iJones et al. (2002). The ratio between free 
spectral range and spectral purity is called finesse and is 
a critical parameter of all TFs. In the case of OSIRIS, 
the finesse is of the order of 50 (see OSIRIS User Man - 
ual), which allows the use of Eq. [2] (I Jones et al.l 120021 ). 
Spectral purity and FWHM coincide when the finesse is 
large. The wavelength at the optical center (Ac in Eq. [IJ 
as well as the FWHM are set by tuning the etalon. In 
our observations the FWHM was set to 18 A (see Table 
1) rendering a free spectral range of the order of 900 A. 

A monochromatic image at a particular wavelength 
(e.g. Ha) over the entire FoV of the detector can be ob- 
tained if we have a sequence of images where Ac between 
successive images increases by <FWHM. For the image 
sequences we have, this criterion was not strictly met (AA 
between successive images was 20 A for a FWHM=18 A). 
The consequence of this is that there are annular zones 



with missing data. However, multiple observations with 
dithered positions helped us to fill-in data for these zones. 
In what follows we describe the method we have adopted 
for reconstructing a monochromatic image from MlOl 
image sequences. 

3.2. Reconstruction Strategy 

The observed count rate F(x, y) at a pixel x, y of an 
image is related to the emitted intensity Ix{x,y) by the 
equation: 



F(x,y) 



J h{x,y)Rx{r)dX 



+ Sky(a;,y) (3) 



where R\{r) (Eq. [2|) is the value of the response curve 
at a distance r = y/{x — xqY + [y — VoY from the 
optical center (a;o,yo), i^ is the conversion factor be- 
tween the intensity and the count rate in units of 
ergcm^^s^^/(counts^^), and Sky{x,y) is the count rate 
from the sky at a pixel x,y. In §6, we describe a proce- 
dure to determine the value of k. Ideally, the sky term 
can be estimated if observations of sky frames are taken 
in each scanned wavelength. Given that this involves 
extra telescope time, such observations were not carried 
out. The dependence of sky on x,y is because of two 
causes (1) an intrinsic spatial variation of the sky, and 
(2) a wavelength dependence of sky emission. The first 
of these can be assumed to be negligible over the FoV 
of OSIRIS and hence the sky variation in the image is 
principally due to the second effect. Given the circu- 
lar symmetry of the variation of wavelength, sky val- 
ues are radially symmetric around the optical center, i.e. 
Sky{x,y) = Sky(r). We determined Sky(r) as the aver- 
age count rate from the sky pixels in annular zones of 
width=FI^i/Af/2. 

In general, I\ includes two contributions: (1) Ivmo, 
emission line flux integrated over the line profile, and (2) 
/cont, the continuum flux density. We first describe the 
method for reconstructing monochromatic images in the 
simplest case of only one emission line entering the filters 
in the entire scanned wavelength range. The emission 
lines from nebular sources are extremely narrow with re- 
spect to the typical bandwidth of the TF, and hence the 
count rate from a line depends only on the value of the 
response function at the wavelength of the line, -Riinc('')- 
Hence the integral can be replaced by a simple multi- 
plication, i.e. J Iuncix,y)R\{r)dX = hne{x,y)Rune{r), 
when line-width << FWHM. On the other hand, 
the term due to the continuum from the source is 
f i^^^,(x,v)R^ir)dx _ Cont(x,2/), where Cont{x,y) is the 
observed count rate in the continuum image. With these 
approximations Eq. [3]can be re- written as: 

F{x^ y) = ^"""^"^' y)R^^^i^) + Cont(a:, y) + Sky(r). (4) 

K 

A trivial manipulation of the above equation gives: 
_ K {F{x, y) ~ Cont(a;, y) - Sky(r)) 

ninc\XTy) — — 1-T (,0J 

R\inc[r) 

By making a simple substitution 

C(x, y) = F{x, y) - Cont(x, y) - Sky(r), (6) 
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where C{x,y) is the sky and continuum subtracted count 
rate at a position x, y of the image, the above equation 
can be re-written as: 



hnc{x,y) 



K. C{x,y) 

Rlinc{r) 



(7) 



Due to the long tail of the response curves, an emission 
line from some regions of the galaxy is registered by more 
than one image of the scan. Hence, in general, simultane- 
ous equations of the kind of Eq.[7]can be written for each 
image. The line contribution in the two adjacent images 
can be combined, with the weights for combining being 
determined by the values of the response functions for 
the line at each pixel. A generalized equation for recov- 
ering the flux from any pixel where the line is registered 
by two consecutive images of a scan: iml and im2, can 
be written as: 

hinc(x,y) = -— -— —7- -— (8j 

[-Rlino (?-)], ™i + [^linc(r-)],„2 

For the sake of compactness, from now onwards R\{r) 
will be denoted simply by Rx- 

3.3. Ha and [S ii] scans 

Due to the long tail of the TF response curve, an emis- 
sion line separated from the target line by ~ FWHM 
contributes non-negligibly to the flux received in the im- 
age. This is the case while mapping the Ha line with a 
tunable filter of FWHM~ 18 A, where the contribution 
of the fianking [N ii] lines to the observed flux cannot 
be neglected. The contribution from an unwanted line 
becomes even more important while trying to recover 
the [N ii] lines, as the Ha line is likely to dominate the 
observed flux, rather than the [N ii] lines. Another inter- 
esting case is that of the [S ll] doublet, where both lines 
contribute in most of the images of a typical TF scan. 
These cases can be easily treated by adding more terms 
— one term such as /iinc(2;,2/)-Riino('') for each emission 
line — to the numerator of the first term in Eq. 21 For 
the scan involving the Ha line, the first term in the nu- 
merator of Eq. U should be replaced by: 



-^lino(a;, y)Rlinc -^ luaRxHa + -^6583^A6583 + /6548^A6548, 

(9) 
where Rxua, Rxesss, and -Ra6548 are the values of the 
response function at a radial distance r from the optical 
center for the observed (not the rest-frame) wavelengths 
of the corresponding lines. 
Eq. [7]for the Ha line can then be re- written as. 



Iiia{x,y) 



KC{x,y) 



RxHa 



^65 



(^A6583 



q-Ra6548) 



(10) 



where we have substituted the value of the intrinsic flux 
ratio of 4^^^ = i. The equation has a term involving 

J65S3 J ^ <= 

-fM3 [j^ ^jjg denominator. Thus in order to reconstruct 
the Ha line image, we need to know a priori the flux 
ratio images of the contaminating [Nii]A6583 line with 
respect to the Ha line. The [N ii] line fluxes, on the other 
hand, require a knowledge of the Ha flux, as can be seen 



by the equation for the [Nii]A6583 hue: 

KC{x,y)- Ilia{x,y)RxHa ._,,, 

I6583{x,y)^ . (11) 

-n'A6583 

In this equation, we have neglected the [N ii] A6548 line 
contribution given that it contributes less than 2% in im- 
age sections tuned to maximize the [N ii]A6583 emission. 
Thus, in order to obtain [Nii]A6583 flux map, we need 
to know the Ha flux. An iterative procedure involving 
Eqs. [To] and [ll] is necessary for an accurate recovery 
of both Ha flux and -j^^ ratio. Given the exploratory 
nature of the present study, we have calculated the Ha 
fluxes by fixing -^^^ = 0.1 for all regions. This approxi- 
mation would result in an overestimation of the Ha fiuxes 
by > 2.5%, and underestimation of the -j^sa ratios by 

~ 0.1 for regions having intrinsic -^P^ > 0.2. 

The reconstruction equations for the A6716 and A6731 
lines from the [S ii] scan are: 



'6716 



(a;,2/) 



[KC{x,y)]^ 



[KC{x,y)]^ 
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(12) 



'6731 
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[KC{x,y% 



[KC{x,y)] 
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T>iral 
^A6731 



R 



+ 



3(x,y) 



A6731 ^ l673i{x,y) 



i^A6716 + R 



A6716) 

(13) 



In these equations, iml and im2 denote the image sec- 
tions tuned to maximize the A6716 and A6731 lines, re- 
spectively. Note that, in order to recover the /67i6 image, 
we need the /673i image and vice versa. Unlike the case of 
Y-^ , no apriori values for the intensity ratios of the [S ii] 
lines could be used, given that this ratio is very sensitive 
to the electron density of the regions. We hence resolved 
Eqs. [12] and [131 by assuming a value of 7^^^, and itera- 
tively changing that value until the value at every pixel 
stabilizes within 10% in two successive iterations. We 
checked that the resulting images are the same irrespec- 
tive of the starting value of 4^ 



-f6731 ' 

4. IMPLEMENTATION OF THE METHOD 

We developed a script under the IRAF^ environment 
to implement the method in a user-friendly way. The 
first step of the reconstruction process is obtaining wave- 
length vs. pixel look-up images for every image of the 
scan. These images are created using Eq.[T]and the astro- 
metric solutions as described in §3.1. Values of Ac^" and 
FWHM are taken from Table 1. For each of the emission 
lines, we then created 2-dimensional response images us- 
ing Eq. [2] The procedure we followed is illustrated using 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement mth 
the National Science Foundation. 

^^ Wavelength calibration lamps were not supplied in this initial 
observing run, and hence we checked/recalibrated the value of Ac 
by comparing the filter-convolved SE)SS spectra of 7 H II regions 
to the profiles of the corresponding H II regions formed using the 
observed fluxes in successive images. Error in Ac using this method 
is found to be < 2A, which is better than what could be achieved 
using the sky rings. 
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Fig. 1. — (Top) An illustration of the variation of the wave- 
length over the OSIRIS FoV as expected from Eq. [l] with Ac = 
6588 A (observed wavelength of [Nll]A6583 line in MIDI) and 
FWHM=18 A. The positions where [Nn]A6583, Ha, [Nn]A6548, 
lines have maximum response are marked. (Bottom) The response 
curves, as calculated using Eq. [2] for each of these lines are shown. 
The part of the curves above the horizontal line {rjn^g = 0.4) can be 
used to reconstruct the monochromatic images in the correspond- 
ing lines. 



1-d cuts on one of these images in Fig. [T] Specifically, 
we chose the image with Ac = 6588 A for illustration, as 
this image is capable of detecting three different emis- 
sion lines at different radial bins. The bottom panel 
shows the expected value of the response functions for 
the three lines as a function of the radial distance from 
the optical center. In cases such as this, where more 
than one emission line is registered within the 8' x 8' 
FoV (r = 4'), a second line starts contributing signifi- 
cantly before the line of interest drops to below ^50% 
response, as illustrated in the bottom panel of Fig. [1] It 
is possible to isolate the contribution of each line by se- 
lecting only those zones where the line of interest has a 
response value above a certain level, as described below. 

4.1. Choosing the value for the response cut-off 

The part of a TF image that can be considered 
monochromatic is decided by the value for the response 
cut-off (r/ii,ie) parameter. By carefully choosing r/ii,io, 



it is possible to reconstruct a monochromatic image in 
the line of interest even in the presence of contami- 
nating lines. Only those pixels with a response value 
.Kline > ''yiinc wiU be Considered good for reconstructing 
the monochromatic image in that line. All these pixels 
belong to an annular zone of particular width (circle in 
the case of a line falling close to the optical center). For 
example, with r^unc = 0.4, the image sections for recon- 
structing the Ha and [N ii]A6583 lines in the Ac — 6588 A 
image correspond to the radial zones where the curve lies 
above the horizontal line in Fig. [1] 
Two guidelines are useful for setting the value of Tyune: 

(1) the response for the contaminating line in the selected 
image section has to be less than that for the main line, 

(2) there are no annular gaps in the final reconstructed 
image. The first condition depends on the wavelength 
difference between the contaminating lines, whereas the 
difference between central wavelengths of successive im- 
ages of the scan (A Ac) determines the fulfillment of the 
second condition. For the case of Ha and [Nii]A6583, 
the response value would be the same at a wavelength 
mid-way between the two lines (i.e. 6573 A at rest-frame 
or 6578 A for our MlOl scan given its recession velocity 
of 214 kms'^). At radial zones where Xr < 6578 A, the 
observed count rates would be predominantly from the 
Ha line, implying r^iino ^ 0.45. On the other hand, the 
second criterion requires that A^ be separated from Aho 
by at least AAc/2, implying Tyunc < 0.45 for our MlOl 
scan. Thus rjunc = 0.45 is the optimal value for our 
dataset. The value of rjn^ic could be marginally higher 
than this if dithered images are able to fill in data-less 
annular zones, or lower if the contaminating line con- 
tribution could be subtracted to within a few percent 
accuracy using an iterative procedure. For example, for 
our MlOl scan, dithering between the images was suffi- 
cient to be able to fill-in the data gaps for rjunc = 0.5. On 
the other hand, the SNR of the [Nii]A6583 images was 
not good enough to iteratively subtract the [Nii]A6583 
contamination in the predominantly Ha pixels. The user 
would be required to select the value of rjn^c based on the 
data parameters and the specific scientific objective. 

4.2. Coadding monochromatic image sections 

In Fig. [21 we show the response for the Ha line in 
consecutive images of a TF scan. For a value of r^iinc = 
0.4, in certain ranges of radial zones, the line is registered 
in two consecutive images. Thus, for ryunc = 0.4 there 
is redundant data for some pixels for the construction 
of the monochromatic images. This redundancy can be 
used to our advantage to get deeper images, by coadding 
the pixel values from both images that contribute to that 
zone. The response curves are also coadded to get a net 
response curve such as shown by the solid line in Fig. [2l 
The coadded line image is divided by the net response 
image to get the entire image in the same flux scale (see 
Eq.©. 

4.3. Preparation of the continuum and sky images 

Fundamental to the reconstruction process is the deter- 
mination of pure count rate in the line C{x,y) from the 
observed count rate F(x,y). This is carried out through 
Eq. [SI which involves the estimation of sky and contin- 
uum images. The sky count rates Sky(r) are calculated 
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tral wavelength depends directly on the first derivative 
of the response function (Eq. [2]) with respect to A. 



1 a 3 4 

radial distance (arcinin) 

Fig. 2. — The net response curve for the reconstruction of the 
Ha image for »7iinc = 0.4 is shown by the sohd curve. Ha response 
curves in individual TF images are shown by the dashed hues. The 
central wavelength Ac (in A) for each image is indicated above the 
corresponding response function. 



as the median values in semi- annular zones, with a width 
equal to FWHM/2 = 9 A. 

The bluest wavelength of our scans is Ac = 6528 A, 
which is the filter with least contamination from any 
emission line. We used this image as a first-guess image 
of the continuum. Even for this filter, the Ha line con- 
tributes more than 3% for pixels closer than 90" from 
the optical center. We defined a parameter Tycont such 
that only those pixels for which the response to detect 
any line does not exceed rycont are considered for contin- 
uum (i.e. only pixels with Rune < ??cont)- For example, 
for the image with Ac = 6528A, with ?7cont = 0.03 only 
pixels that are more distant than 90" from the optical 
center satisfy this condition. Even for these pixels, we 
estimated the line contribution using the reconstructed 
emission line maps, and subtracted it from the continuum 
pixels. Individual sections contributing to the continuum 
are stitched together (averaged when more than one im- 
age contributes to a pixel) to obtain the final continuum 
image. This image is input as the new continuum image 
and the entire process is repeated. From the simulated 
data (see §5.3, Run 1), we found that the continuum val- 
ues converge after 3 iterations. A small value of ?7cont 
would result in various annular zones where there are no 
pixels satisfying the condition i?iino < f?cont- From sim- 
ulated data, we found that line-free continuum images 
could be obtained even for Tycont as large as 0.2. 

An IRAF package containing the scripts developed as 
part of this work is available to users upon request and 
will be available for downloading fr om the LUS webpage 



http://www.inaoep.mx/~gtc-lus/ 



5. FLUX ERRORS IN THE RECONSTRUCTED IMAGES 

5.1. Flux error due to tuning error 

The central wavelength of the OSIRIS TF can be set 
only to an accuracy of 1 A (See OSIRIS TF User Man- 
ual) . This uncertainty in tuning the TF leads to an error 
in the recovered emission line fluxes. The error {SF) in 
flux (F) due to an uncertainty of (5 Ac in setting the cen- 



SF 



SR 



1 fdRy 



Substituting the value of ^r^, we get 



dRx 



SF 

F 



= 4Rx 



FWHM j 



R\ 



-1 



(15) 



Maximum error is introduced for image pixels where the 
TF response R\ for the line of interest is 0.5. Thus, flux 
errors in the reconstructed image can be reduced if we 
use ryiine > 0.5. It may be recalled from the analysis in 
§4.1 that the lower limit of Tyiino depends on A Ac. In 
order to have Tynno > 0.5, AAc should be less than the 
FWHM. 

We used the Eq. [15] to calculate the errors in the line 
fluxes, and summed the errors in quadrature to calcu- 
late the errors in the flux ratios, for nominal value of 
(5Ac = 1 A. Our results are summarized in Fig. [S] In 
the three panels, we plot the errors in the calculated 
Ha fluxes (top), the [Nii]A6583/Ha ratio (middle) and 
[S ii]A6717/6731 ratio (bottom) as a function of the sam- 
pling parameter, for two extreme values of FWHM per- 
mitted by the TF. The value of ryuno is calculated in such 
a way that there are neither gaps nor overlapping pixels 
in the reconstructed image for a given sampling AAc. 

The most notable characteristic in these plots is that 
the errors in all these three quantities are lower for larger 
values of FWHM. This may seem counter-intuitive, and 
is the result of a sampling error of 1 A being a larger 
fraction of FWHM = 12 A than that for FWHM = 
18 A. In other words, the TF response function falls less 
steeply for larger FWHM, thus resulting in little errors 
as compared to that for smaller FWHMs. As expected, 
finer sampling results in smaller errors on the derived 
quantities, especially when the sampling rate is less than 
0.7 X FWHM. 

Line fluxes can be calculated with better accuracy than 
the flux ratios. The error in the [Nii]A6583/Ha flux 
ratio takes into account the contamination by Ha in the 
[Nii]A6583 filters. This cross talk also makes the errors 
on the [Nii]A6583 line flux marginally larger than those 
for the Ha flux. The calculated errors on the individual 
[S ii] line fluxes are similar to that for the Ha flux. 

It may be noted that a recessional velocity of 45 km s~^ 
produces a shift of lA in the wavelength. Thus, in re- 
gions where kinematic deviations of this order or more 
are likely to exist (e.g. nuclear regions of galaxies), kine- 
matical data are required for an accurate flux calibration. 
In the absence of such data, eq. 15 can be used, where 
(5 Ac is to be replaced by the uncertainty in the wavelength 
due to Dopper effect, to estimate the flux errors due to 
kinematic deviations. 

5.2. Reconstruction of simulated TF images 

The dataset we have for MlOl is reconstructed with 
'7iino = 0.5, implying maximum errors of around 11% in 
the recovered Ha fluxes according to Fig. [S] In order to 
check the implementation of the equations of §3 in our re- 
construction scripts, we need datasets with much smaller 
errors. We hence created artificial datasets by simulating 
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Fig. 3. — Errors in flux and flux ratios due to an uncertainty of 
(5 Ac = 1 A in setting the central wavelength of TF, as a function 
of user-selected TF parameters A Ac and FWHM. The errors are 
expressed as percentage values of the plotted quantity, (top) Errors 
on the Ho flux, (middle) errors on the [Nll]A6583/Ha ratio, and 
(bottom) errors on the [S ll]A6717/6731 ratio. The scale on the 
top gives the maximum value of response cut-off (»?iine) that the 
data permit for given values of AAc and FWHM, without having 
data gaps in the reconstructed image (see §4.1 for details). Errors 
expected in our dataset of MlOl for an assumed error of <5Ac = 1 A 
are marked by solid circles. 



the observations of an extended emission line source with 
a scanning TF, making use of Eqs. 1 and 2. The simu- 
lated dataset also allowed us to quantify the accuracy of 
reconstruction to various instrumental parameters that 
have non-zero uncertainties. Two scans were simulated: 
(1) TF with a FWHM=18 A, Ac = 6528 + (i - 1) x 20 A, 
where i varied from 1 to 9, (2) TF with a FWHM=12 A, 
Ac = 6528 -I- (i — 1) X 10 A, where i varied from 1 to 
18. The reddest wavelength in both cases is 6688 A. The 
FWHM in the two settings are among the extreme val- 
ues permitted with the red etalon, with the former one 
simulating the dataset we have for MlOl. The intensity 
of the extended source is set constant over the FoV of 
OSIRIS with I(Ha)=l. The [Nii]A6583 line is fixed at 
10% of that of Ha, which is amo ng the range of val- 
ues observed in Hii regions (Denic olo et al.ll2002[ ). The 
[Nii]A6548 line intensity is fixed at one third of that of 
the [Nii]A6583 line. We also added continuum sources 
at some fixed positions on the image, in order to test 
the accuracy of the continuum subtraction. The intensi- 
ties of the continuum sources are adjusted such that the 
Ha emission equivalent widths at the positions of the 
continuum sources are 10 for some sources and 100 for 
others. Gaussian noise of two different values is added 



to the images, resulting in two sets of simulations: the 
high and low Signal-to-Noise Ratio (SNR) images cor- 
respond to SNRs of 100 and 3 for the [Nii]A6583 line 
image, respectively. A line-free continuum image is also 
generated in the simulated set. The lines were redshifted 
by 214 kms~^, corresponding to the recession velocity of 
MlOl. 

5.3. Parameters that limit the reconstruction accuracy 

Reconstruction with ideal datasets (Run 0): The equa- 
tions described in §3 should allow us to recover the fiuxes 
of the Ha and [Nii]A6583 lines, in the absence of un- 
certainties in various parameters that control the TF 
imaging. We simulated this case by using high SNR 
data and subtracting the simulated line-free continuum. 
The errors in recovering the Ha fiux and the fiux ratio 
[Nii]A6583/Ha are as smaU as 0.03% and 0.2% respec- 
tively. This test establishes that the equations developed 
in §3 are correctly implemented in the script. 

However, there are many sources of error in real ob- 
servational datasets, especially when trying to maximize 
the available observing time. We investigate the accu- 
racy in the recovered flux due to the following 3 sources 
of error: (1) the continuum images are not completely 
free of emission lines, (2) the sky images are obtained 
by in-frame object-free pixels, not from an off-source sky 
image, (3) different images within a scan have non-zero 
ditherings. We discuss each of these cases in detail below. 

Line contamination in continuum images (Run 1): 
Due to the long tail of the TF response, there is non-zero 
line contribution in all TF images, including those cen- 
tered bluer than the bluest line in a scan (e.g., Ac = 6528 
is only - 1 FWHM blueward of the [N ii]A6548 fine, and 
< 2 FWHMs of the bright Ha line). In order to obtain 
line-free continuum fluxes from the sequence of images 
we have, we followed the iterative procedure described in 
§4.3. We obtained a continuum image from the sequence 
of simulated images, and compared it to the simulated 
pure continuum image. The difference between the two 
images is found to be as small as 0.1% after 3 iterations. 
Hence, errors on the fluxes in our reconstructed images of 
MlOl are not dominated by the limitations of obtaining 
line- free continuum images. 

Non-uniform sky (Run 2): Errors in sky value subtrac- 
tion and flat-field corrections introduce residuals in the 
sky value locally. We parametrized this error in terms of 
the rms noise value (cr) of this image. The residual sky 
value in different parts of the image is allowed to vary 
between —la and -|-lcr, for the simulated set of images. 
The [Nii]A6583/Ha ratios are affected by 5% in pixels 
where [N ii]A6583 is detected with a SNR of 3. Thus, for 
pixels with SNR> 3 for the fine of [N ii]A6583, sky sub- 
traction is not a serious problem in our images of MlOl. 

Consequences of Image-dithering (Run 3): It is a nor- 
mal practice to dither images between any two repeat 
observations in order to avoid detector blemishes spoil- 
ing any interesting feature. These dithered images are 
registered to a common coordinate system before com- 
bining them. In a TF observation, the optical centers 
in the dithered images correspond to different astromet- 
ric coordinates and hence the corresponding pixels in the 
registered images do not have the same wavelength (see 
Eq. 1). Thus, in the combined image, a given pixel has 
contributions from marginally different wavelengths. We 
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parametrized this effect as an error in the optical center 
Src- We studied the reconstruction accuracy for vari- 
ous values of Sr^ between 1" and 5". A dithering of 1" 
between different images of a scan can produce errors 
of the order of ^ 5% in the recovered Ha flux for the 
FWHM=18 A fihers. The recovered [Nii]A6583/Ha ra- 
tio lies between 0.05-0.15 (50% error over the simulated 
value of 0.1) for this case. Combining images that are 
dithered by more than 1" would introduce more than 
20% error on the recovered Ha flux, making it unusable 
for most applications. In the case of our observations of 
MlOl, the dithering within a scan was less than 0.25", 
and hence the error in the recovered Ha flux due to as- 
trometric registration of all images of a scan (say PI) is 
less than a few percent. 

5.4. Recommendations for observing extended sources 

After studying the effect of various parameters on the 
fluxes in the reconstructed images, we find that the 
maximum error on our dataset for MlOl arises due to 
the present uncertainty of ~ 1 A in setting the central 
wavelength of a TF observation. This uncertainty af- 
fects more the images taken with narrower TF observa- 
tions, with the accuracy of Ha fluxes being ~10% and 
of [Nii]A6583/Ha ratios of - 18% for the TF images 
with FWHM = 18 A, as can be inferred from Fig. El 
Corresponding errors with FWHM = 12 A are 16% and 
25%, respectively. Errors may be reduced by carrying out 
the scan with a finer saniphng. With FWHM = 12 A 
a sampling of around 5 A would be required to reduce 
the error levels to ^ 10%, that can be achieved with 
15 A sampling with FWHM = 18 A. Thus a factor of 
3 more exposure time would be required to map a given 
emission line over the entire FoV of OSIRIS with the 
narrower TF. However, for low SNR pixels or regions 
for which sky and/or continuum subtraction errors con- 
tribute more than the plotted errors {SNR < 10), the 
errors are expected to increase proportionally with the 
FWHM, and hence observations with FWHM = 12 A 
would have an advantage by a factor of 1.5 over those 
with FWHM — 18 A. Thus, if the interest is in detecting 
d iffuse faint emissions, a F WHM = 12 A is preferable. 

iLara-Lopez et al.l (|2010( ) carried out simulations to de- 
termine the best combination of FWHM and sampling 
(A Ac) for optimal emission line flux determinations of 
emission-line galaxies of redshifts between 0.2 and 0.4 
with GTC/OSIRIS. They found that a FWHM of 12 A 
and a sampling of 5 A are the optimal combination that 
allows deblending Ha from the [Nii]A6583 line with a 
flux error lower than 20%. It is relevant to note that in 
their simulations, the fluxes were not corrected for the 
response curve of the TF, and the quoted flux errors are 
due to the unavailability of the redshifts of the detected 
galaxies. Hence, it is natural to expect lesser errors for 
the narrower fllters. 

The next source of error comes from image dithering. 
For an image dithering of 1", the errors could be as large 
as 5% for Ha and 50% for [Nii]A6583/Ha. For larger 
values of dithering, errors on [N ii]A6583/Ha are unrea- 
sonably high. It is advisable to have the same telescope 
position (dithering < 1") for all the images constituting 
a given TF scan. If more than one scan is available for a 
fleld with dithering of more than 1" between the scans. 



it is advisable to reconstruct an emission-line image from 
each scan, and then register and combine them. We rec- 
ommend that a TF scan intended to obtain a monochro- 
matic image in an emission line should have one image 
observed at least one FWHM blueward of the bluest 
line in the TF sequence, to facilitate accurate continuum 
subtraction. 

6. FLUX CALIBRATION OF MONOCHROMATIC IMAGES 

We have explored a new procedure to flux calibrate 
the TF images using the in-frame fleld stars. The proce- 
dure, in principle, can be applied to any optical narrow- 
band image of a fleld that contains stars with photom- 
etry in the optical broad-bands. In particular, we used 
the SDSS photometry in griz bands of stars in the field 
of MlOl, and the recently released SDSS Stellar Spec- 
tral d a,tabase of stars coverin g a wide range of spectral 
types ()Abazaiian et al.l l2009l. These spectra are of me- 
dian resolution in the 3500-9000 A wavelength range. In 
order to have a good spectral library we selected, from 
the original Stellar Spectral database, only those stars 
with g =14-18 mag having no gaps, jumps or emission 
lines in their spectra. 

The general procedure involves obtaining a stellar spec- 
trum that fits the griz magnitudes of the stars in the 
field of the TF image. Basically, we are using the stel- 
lar spectrum to interpolate the broad-band fluxes at the 
wavelength of the TF. The spectra are integrated in the 
gri z-haiids to get their synthetic magnitudes which are 
then fitted to the griz magnitudes of a star in the field 
of interest; both Spectral Energy Distributions (SEDs) 
are normalized at the r-band. 

The best-fit spectrum is chosen by minimizing the x^ 
obtained in the 4 bands. The spectrum of the best-fit 
SDSS star is then de-normalized by multiplying it by the 
r-band flux of the field star. The resulting spectrum is 
convolved with the response function of the OSIRIS TFs 
(Eq. [5]) to obtain a smoothed spectrum of the field star. 
The flux at A^, where r is the radial distance of the star in 
the OSIRIS fleld, is multiplied by the effective bandwidth 
of the TF, to estimate the flux intercepted by the TF. 

We then carried out aperture photometry of the se- 
lected SDSS stars on all the TF images. This photometry 
is used to obtain observed count rate of each star in each 
TF image. The estimated flux is divided by the observed 
count rate of the star in that TF to obtain the calibra- 
tion coefficient k. The procedure is repeated for all the 
good SDSS stars in the observed field to obtain a set of k 
values. We note that the observed count rate is initially 
corrected for the effects of extinction and the efficiency 
of the order sorter filter, and hence the n obtained from 
different stars in different TFs can be directly compared 
with each other. 

6.1. Relative errors in the Calibration coefficients 

The availability of thousands of stellar spectra covering 
the entire range of spectral types to fit the photometric 
data of fleld stars ensures that there is at least one spec- 
trum that truly represents the spectrum of the field star. 
Generally, there are 10-15 spectra whose griz-hand x^ 
is within 10% of the best-fit spectrum. We obtained a 
mean and rms of these spectra at every sampled wave- 
length. The rms error was found to be less than 1% for 
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Fig. 4. — An illustration of the calibration procedure adopted, where we fit the SDSS griz photometry of in-frame field stars with the 
spectra of stars in the SDSS spectral catalog, (left) A section of the OSIRIS continuum image of MIDI. Four stars with SDSS photometry 
are identified, (right) SEDs in griz bands of the 4 identified stars (asterisk), superposed on the average of 12 best fitting spectra from 
SDSS spectral catalog (line). The rms deviation from the mean of 12 best-fit spectra at each wavelength is shown by the gray band around 
the mean spectrum. The 12 spectra differ by less than 1% (< rms >r denotes the average rms value within the r-band), which makes this 
method very attractive for the calibration of OSIRIS TF data. 



wavelengths between the g and z bands. Thus the rel- 
ative flux calibration for different TF settings is better 
than 1%. The absolute flux error depends on the error 
in the SDSS r-band magnitude of the field stars. An 
illustration of the method followed is shown in Fig. |4l 

In Fig. [5l we show the variation of the calibration co- 
efficient along an Ha TF scan containing 9 images. In 
this figure, two kinds of systematic variations are seen: 
(1) the variation of the coefficient from one TF observa- 
tion to the next, and (2) the variation of the coefficient 
between the two CCDs. 

Variation of calibration coefficient during a TF scan: 
given that the airmass variation from one observation 
to the next is taken into account in obtaining the k, this 
variation was not expected for a photometric night. The 
rms dispersion for the 4 stars that are used to calculate 
K is around 3%, which is much smaller than the overall 
variation (9% over the mean value). The trend of the 
variation is identical for the two scans. 
Variation of calibration coefficient between the two 
CCDs: OSIRIS uses two CCDs, along with 2 separate 
electronics to cover the total field of view. The differ- 
ences in the background levels between the two CCDs 
are easily noticeable. From our analysis, we find that 
there is a 6 ± 3% difference between the efficiencies of 
the two CCDs for the Ha scan. The most likely rea- 




FlG. 5. — Variation of the calibration coefficient k for 9 images 
of the two Ha TF scans of MlOl. We used 8 stars, 4 on each CCD 
covering the complete field of view. The coefficients for the right 
(dashed line) and left (dotted line) CCDs are separately shown. 
The solid line shows the variation without distinguishing on which 
CCD stars are located. In the top panel, we show the coefficients 
for scan 1 (PI) and in the bottom, for scan 2 (P2). 



son for this difference is the value of the gain parameter, 
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Table 2 
Derived values of calibration coefficients 



TF Scan k 


Units 


Ha+[Nii]-Pl 6.54 ±0.27 10" 
Ha+[Nn]-P2 6.63 ± 0.28 10" 
Sii-Pl 7.06 ±0.14 10- 
SII-P2 7.06 ±0.21 10- 
SII-P3 7.27 ±0.23 10- 


-18 ergs-icm-^/(ADU/s) 
-18 ergs-lcm-2/(ADU/s) 
-18 ergs-lcm-2/(ADU/s) 
-18 ergs-icm-2/(ADU/s) 
-18 ergs-icm-2/(ADU/s) 



which is 0.95 electrons/ADU for the right CCD and 0.91 
electrons/ADU for the left CCD (private communication 
from GTC technical staff) during the Ha scan. 

When these systematic variations in k are taken into 
account, the calibration coefficients obtained using differ- 
ent stars in different CCDs and in different TF images, 
agree to within 3% of each other. The mean cahbration 
coefficients obtained from the stars on the right CCD of 
the first image of the Ha and [S ii] scans are given in Ta- 
bloid] The difference in the calibration coefficients for the 
two Ha scans is of the order of 1%, whereas for the [S ii] 
scans this difference is around 3%. The mean difference 
in the cahbration coefficients for the Ha and [S ii] scans 
is '^ 7%. The k factors in the reconstruction equations 
take into account both these systematic variations, thus 
aUowing us to combine the monochromatic images from 
different images of a scan. 

7. COMPARISON OF RECONSTRUCTED IMAGES 

In Fig. m we show the reconstructed Ha emission-line 
image of the south-west part of MlOl. The Ha line be- 
ing a good tracer of ionized gas, this image shows the 
location of H ii complexes and ionized filaments. It can 
be seen from the image that we are able to detect H Ii 
regions in the entire field of view. The 5cr surface bright- 
ness limit in this image is 8x 10""'^'' ergcm"^ s~^ arcsec"^. 
All the regions seen in the image are recovered in the two 
independent scans. In the following paragraphs, we com- 
pare the fluxes of recovered emission lines of H ii regions 
in these two independent scans. Our Ha fluxes are also 
compared with the Ha fluxes obtained using traditional 
narrow-band imaging technique. 

We performed photometry of 23 H ii regions using 
an aperture radius of 4" on our reconstructed Ha, 
[N ii]A6583, [S Ii] images and also on an Ha+[N ii] image 
that was obtained in the traditional way using narrow- 
band filters (downloaded from NED: Telescope: KPNO 
Schmidt; Observers: B. Greenawalt, and R. Walterbos). 
The selected regions are among the brightest regions in 
the FoV. 

With a sampling of A Ac = 20 A for our dataset, r/unc 
should be at least 0.45 in order to obtain an image with- 
out data gaps. However, we found systematically larger 
dispersion in fluxes for regions with response values less 
than 0.5. Therefore we carried out the reconstruction 
with ryiine = 0.5. As a result, eight regions have data 
in only one of the dithered scans. In Fig. [7^, we show 
a comparison of the fluxes of the 15 H ii regions whose 
fluxes could be measured on the reconstructed images 
from both the scans. The mean ratio of the Ha fluxes is 
unity over 2 orders of magnitude in fiux, confirming the 
accuracy of the method we adopted for the fiux calibra- 
tion. The rms dispersion of the Ha fiuxes of the same 



region obtained on images for the two scans is 16%. This 
value is within the flux errors expected for the 1-2 A er- 
ror in Ac- 

In order to check whether the reconstruction process 
rightly reproduces relative fiuxes of different regions over 
the entire FoV, we plot the ratio of the Ha+[N ii]A6583 
fluxes measured from our reconstructed images (PI and 
P2), and that from traditional narrow-band filters in 
Fig. [Tb. Our fluxes are given as the average when data 
are available from both scans. As the intention here is 
a comparison of relative fluxes, we set the mean value 
of the flux ratio to unity. There is ~ 11% scatter on 
this mean value, which is marginally better than that 
between PI and P2, as expected due to the use of av- 
eraged fluxes. There is a marginal trend for the mean 
ratio to be ^ 5% different between the bright and faint 
regions, which is nevertheless smaller than the scatter. 
We did not find any trend of these two ratios against the 
distance of the region from the optical center. 

Three sources of error are included in calculating the 
sizes of the error bars plotted in Fig. [T] They are (1) 
photon noise of the object, (2) the error in the subtrac- 
tion of the sky value and (3) the error in the recovered 
flux due to an error in Ac- The last of these errors, which 
is discussed in detail in §5, dominates for the 23 regions 
for which we performed photometry, contributing around 
10% for the majority of the regions. 

As a second test of the reliability of the fiux ratios ob- 
tained using the TF images, we compare the flux ratios 
of emission lines relevant for diagnostic diagrams with 
those values obtained using spectra from the literature 
for the same regions. Seventeen H ii regions in MlOl were 
obser ved spectroscopically by SDSS (jSholudchenko et all 
l2007t ) . Six of these lie within the usable FoV of our im- 
ages (radial distance from the optical center < 3.75'). 
Flux ratio of nebular diagnostic lines from our recon- 
structed images are compared with those obtained from 
the SDSS spectra in Fig. [51 Our fiuxes were obtained 
over apertures of 4" radius at the coordinates associated 
with the SDSS spectra. The adopted apertures, though 
almost 3 times bigger than the fiber sizes of SDSS spec- 
tra (3" diameter), are the minimum area over which re- 
liable fiuxes can be measured in our images. Given that 
the spectroscopic ratios of giant H li regions are not ex- 
pected to vary much with aperture size, our relatively 
bigger apertures are not expected to introduce additional 
errors. The difference in the fiux ratios between ours and 
SDSS values are plotted both against flux ratios (in the 
left panel) and Ha fluxes (in the right panel). The plot- 
ted error bars take into account all the errors discussed 
in the paragraph above. The errors in the spectroscopic 
ratios are expected to be almost negligible, and hence we 
did not include these errors in our analysis. The major- 
ity of the points lie close to the horizontal line within the 
plotted errors, indicating that the ratios of lines relevant 
for diagnostic purposes can be obtained using TF imag- 
ing. There are a few ratios that deviate from the spec- 
troscopic ratios by more than the estimated errors. The 
most important of them is the right-most point in panel 
(a). Apart from being the faintest in Ha, the spectro- 
scopic [N ii]A6583/Ha value for this region is 0.4, a value 
that is too high as compared to the assumed value of 0.1 
in Eqn [TO] This is the most likely reason for the large 
deviation of this region. Kinematics of individual regions 
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Fig. 6. — Reconstructed Ha image of the south-west part of MlOl. Several Hll regions, including the massive Hll complex NGC 5447 
(a = 14 : 02 : 30, 5 = 54 : 16 : 15), as well as several filamentary structures can be seen in this image. The brightness of the Ha structures 
is shown by the gray bar at the bottom of the image in units of ergs~^ cm~^ arcsec"^, spanning a range from —Sir to 512(T in logarithmic 
scale. 



with radial velocities > 25 kms~^ (i.e. 0.5 A around the 
Ha wavelength) can also be responsible for the observed 
deviations of some of the points. The errors in the mea- 
sured radial velocities of the regions from SDSS spectra 
do not permit us to carry out a more elaborate analysis 
of the deviations of individual regions. 

In summary, the ratios recovered from the TF images 
for individual regions have larger errors than the corre- 
sponding spectroscopic ratios. Nevertheless, the capabil- 
ity of TF images to derive such ratios for large number 
of regions, makes TF imaging a scientifically attractive 
option, as we illustrate below using the flux ratio of 23 
H II regions in MlOl. 



In Fig. H we plot the [N ii]A6583/Ha, [Sii]A6716 -I- 
A6731/Ha, and the [S ll]A6716/A6731 ratio of Hii re- 
gions against their Ha fluxes in panels (a), (b) and (c), 
respectively. In each of these panels, we indicate the ob- 
served/theoretical ranges of these ratios in H ii regions, 
which cover ~ 10 times the errors on the ratios. This 
relatively large dynamic range, combined with the in- 
trinsic multi-object capability of TF imaging, makes it 
competitive with traditional spectroscopic observations. 

The [N ii] and [S ii] lines originate in low ionization 
zones in H ii regions. Their ratio to the Ha flux depends 
not only on the ionization parameter, but also on the 
abundance of these ions. The utility of these ratios for 
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Fig. 7. — Relative errors on the Ha fluxes of 23 H II regions in 
MlOl. (a) Ratio of the Hq fluxes measured on the reconstructed 
image for the dithered position PI to those for position P2; (b) 
ratio of the H II region fluxes measured on an Ho+[Nll] image 
taken from NED to those on our images, both plotted against our 
Ha fluxes. The dotted horizontal lines denote Icr scatter over the 
mean ratio. 
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Fig. 8. — Comparison of important diagnostic ratios from 
our reconstructed images with those obtained using SDSS spec- 
tra for 6 regions in common. In the left panels the difference 
between the ratios (e.g. A([Nll]/Ha) = ([N II]A6583/HQ)^^J.^- 
([NII]A6583/Ha)J,JJ,J,) is plotted against the SDSS ratios, whereas 
in the right panels, the difference between the ratios is plotted 
against our Ha fluxes. In general, the diagnostic ratios are repro- 
duced within the ranges allowed by the estimated errors. See text 
for more details. 
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Fig. 9. — Flux ratios of nebular diagnostic lines plotted against 
the Ha fluxes for the 23 H II regions of Fig. [T] In panels (a) and 
(b), we show [Nll]A6583/Ha and [S ll]A6716 -f A6731/Ha, whereas 
in panel (c) we show the density-sensitive ratio [S ll]A6716/A6731. 
The upper limits o bserved in extragaiactic H II regions from 
IDenicolo et al.l II2002I ) are indicated by the horizontal dotted lines 
in the top two panels, whereas t he theoretically valid range for 
the [Sll] line ratios from [Osterbro ck & Feiiand (2006) is shown in 
panel (c). All the regions have observed ratios in the range ex- 
pected for H II regions, illustrating the capability of TF imaging 
for obtaining these diagnostic ratios. 



diagno stic purposes has been discussed in [Baldwin et alJ 
(1981). In Fig. [ini we show the values for these ra- 
tios expected for a range of ionization pa rameters (R) 
and m etallicities (Z), using the models of iDopita et alJ 
()200(3[) . The observed values of these ratios in our se- 
lected 23 regions lie within the range of model val- 
ues. The regions having high values of [Nil] /Ha have 
most likely hig h nitro gen abundances, as was found by 
IDenicolo et al.l (J2002D for H Ii regions in a sample of 
nearby galaxies. Thus, OSIRIS TF imaging is very 
promising for the study of nebular line ratio diagnostics 
of nearby large galaxies. 

8. SUMMARY 

In this work, we have explored the capability of tun- 
able filter imaging with the OSIRIS instrument at the 
10.4-m GTC, using real data. The changing wavelength 
across the field and the non-flat functional form of the 
response curves makes it essential to have a sophisticated 
analysis package to completely take advantage of the rel- 
atively large field of view of the instrument. With the 
set-up that we have used to observe MlOl, we were able 
to obtain monochromatic images in the emission lines of 
Ha, [Nii]A6583 and the [S Ii] doublet. We demonstrate 
that line fiuxes and their ratios for prominent H II regions 
can be obtained to better than ^15%, which is basically 
limited by the current 1 A uncertainty in setting the cen- 
tral wavelength of the TF. Though these errors are much 
larger than the spectroscopic ones, the multi-object ca- 
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Fig. 10. — Nebular diagnostic diagram involving [Nil] A6583/Ha 
and [Sll]A6716/A6731 for Hll regions in MlOl. All the observed 
points lie within locii of models of D opita et al., ^2006 ) defined by 
metallicities between OAZQ and 2Z0, and the logarithmic ion- 
ization parameter log R between 2 (low pressure regions) and —6 
(high pressure regions). 



pability of TF imaging, combined with the relatively high 
sensitivity of the flux ratios to variations in density, ex- 
citation and metallicity, makes TF imaging an attractive 
option for investigating the point-to-point variation of 
these physical quantities in galaxies. We also demon- 
strate that the emission-line maps can be flux-calibrated 
to better than 3% accuracy, using the griz SDSS mag- 
nitudes of in-frame stars, and the spectral database of 
SDSS, without the need to invest extra telescope time to 
perform the photometry of spectrophotometric stars. 
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